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ABSTRACT 



The workshop on "Hadron-Hadron and Cosmic-Ray Interactions at multi-TeV Energies" held at the 
ECT* centre (Trento) in Nov.-Dec. 2010 gathered together both theorists and experimentalists to 
discuss issues of the physics of high-energy hadronic interactions of common interest for the par- 
ticle, nuclear and cosmic-ray communities. QCD results from collider experiments - mostly from 
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the LHC but also from the Tevatron, RHIC and HERA - were discussed and compared to various 
hadronic Monte Carlo generators, aiming at an improvement of our theoretical understanding of soft, 
semi-hard and hard parton dynamics. The latest cosmic-ray results from various ground-based obser- 
vatories were also presented with an emphasis on the phenomenological modeling of the first hadronic 
interactions of the extended air-showers generated in the Earth atmosphere. These mini-proceedings 
consist of an introduction and short summaries of the talks presented at the meeting. 
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Introduction 

The origin and nature of cosmic rays (CRs) with energies between 10'^ eV and the so-called Greisen- 
Zatsepin-Kuzmin (GZK) cut-off at about 10^° eV HI, recently measured by the HiRes [2] and Auger 131 
experiments, remains a central open question in high-energy astrophysics with very interesting con- 
nections to particle physics and, in particular, to Quantum-Chromo-Dynamics (QCD) at the highest 
energies ever studied. One key to solving this question is the determination of the elemental compo- 
sition of cosmic rays in this energy range. The candidate particles, ranging from protons to nuclei 
as massive as iron, generate "extensive air-showers" (EAS) in interactions with air nuclei when en- 
tering the Earth's atmosphere. The determination of the primary energy and mass relies on hadronic 
Monte Carlo (MC) models which describe the interactions of the primary cosmic -ray in the upper 
atmosphere. 

The bulk of particle production in such high-energy hadronic collisions can still not be calculated 
within first-principles QCD and general principles such as unitarity and analyticity (as implemented 
in Regge-Gribov theory) are often combined with perturbative QCD predictions for high-/7r pro- 
cesses, constrained by the existing collider data {Ejah ^ 10^^ eV). Important theoretical issues at 
these energies are the understanding of diffractive and elastic hadronic scattering contributions, the 
description of hadronic forward fragmentation and multi-parton interactions ("underlying event"), 
and the effect of high parton densities ("gluon saturation") at small values of parton fractional mo- 
mentum X = PpamnI Pproton- Indccd, at thcsc energies, the relevant Bjorken-x values are as low as 10^', 
where effects like gluon saturation and multi-parton interactions, particularly enhanced with nuclear 
targets, are expected to dominate the early hadron collision dynamics. 

The current energy frontier for hadron collisions in the laboratory is reached at the Large Hadron 
Collider (LHC), currently under operation at CERN. The measurement of inclusive hadron produc- 
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tion observables in proton-proton, proton-nucleus, and nucleus-nucleus collisions, at LHC energies 
(equivalent to Eiat ~ 10^' eV) will provide very valuable information on high-energy multiparticle 
production, and allow for more reliable determinations of the CR energy and composition around 
the GZK cutoff. In the high luminosity phase of LHC, each bunch crossing will lead to several 
proton-proton interactions, increasing even more the importance of understanding the background 
from diffractive and soft particle production. Semi-hard particle physics will allow one to test the 
boundaries of the applicability of perturbative QCD in the region where low-x gluon saturation phe- 
nomena become increasingly important and may even dominate particle production. 

All LHC experiments feature detection capabilities with a wide phase-space coverage without paral- 
lel, in particular in the forward direction, compared to previous colliders |4|. Such capabilities allow 
for a (fast) measurement of global hadron-hadron collision properties (inelastic - including diffrac- 
tive - cross sections, particle multiplicity and energy flows as a function of pr and pseudorapidity, ...) 
even with the moderate statistics of a first pp and PbPb run. 

The aim of the Workshop was to discuss theoretical and experimental issues connected to hadronic 
interactions of common interest for high-energy particle and cosmic-ray physics. With the recent 
high-quality cosmic-ray results from the HiRes and Auger experiments and the first available LHC 
data it seemed a timely moment to have such a meeting in autumn 2010. The Workshop brought 
together experts, both theorists and experimentalists, in QCD and cosmic-ray physics in view of 
expanding the mutually beneficial interface between two communities currently exploring the physics 
of strong interactions at the highest energies accessible. The talks and discussions on various topics: 

• QCD predictions for high-energy multiparticle production and their implementation in hadronic 
Monte Carlo generators: PYTHIA, HERWIG, SHERPA, PHOJET, DPMJET, QGSJET, SIBYLL, 
EPOS, QGSM, FLUKA; 

• theoretical and experimental developments on diffractive and elastic scattering at high energies; 

• theoretical approaches of multi-parton dynamics and underlying event in hadronic collisions; 

• theoretical and experimental developments on low-;c QCD and inclusive particle production; 

• theoretical developments on modeling of cosmic -ray showers; 

• latest experimental QCD results at colliders: LHC (ATLAS, CMS, ALICE, LHCb, TOTEM, 
LHCf), Tevatron, RHIC, and HERA; 

• latest cosmic-rays measurements in the lO'^-lO^^'eV range: Auger, HiRes, TA, Kascade- 
Grande, Argo, IceCube; 

were organized around four main blocks: 

1. Hadronic collisions at multi-TeV energies: Experimental results 

2. Hadronic collisions at multi-TeV energies: Theory 

3. Cosmic -rays at Ultra-High Energies: Experimental results 

4. Cosmic -rays at Ultra-High Energies: Theory 
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These mini-proceedings include a short summary of each talk including relevant references, the list 
of participants and the workshop programme. We felt that such a format was more appropriate than 
full-fledged proceedings. Most results are or will soon be published and available on arXiv. Most of 
the talks can also be downloaded from the workshop website: 

http://www.cern.ch/CRLHC10/ 

Support of the European Community Research Infrastructure Action under the FP7 'Capacities Spe- 
cific Programme' is acknowledged. We thank the ECT* management and secretariat, in particular 
Cristina Costa, for the helpful cooperation prior and during the workshop and all participants for their 
valuable contributions. We believe that this was only the first workshop of this kind and look forward 
to similar meetings in the future. 



David d'Enterria, Ralph Engel, Torbjorn Sjostrand 
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1 Hadronic collisions at multi-TeV energies: Experiments 



First QCD results from the ATLAS Collaboration 

Mark Suttoi^( University of Sheffield) 

Since the first LHC operation in Nov. 2009 the ATLAS experiment has collected data at sev- 
eral proton-proton centre-of-mass energies with an integrated luminosity of 45 pb^ ' in collisions at 
7 TeV. This large sample has enabled many QCD-related analyses to be performed spanning the en- 
tire kinematic range from soft QCD in minimum bias interactions [5,6| through to the study of the 
underlying event [7] in events with at least one energetic track, through to production of jets with high 
transverse momentum pj lEllll. Events with a prompt photon or a Z or IV have also been observed 
copiously Unidll as have events where the gauge boson is produced in conjunction with a high-pj- 
jetm. 
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Figure 1 : a) The particle azimuthal separation with respect to the leading particle in MB interactions, 
b) The inclusive single jet cross section, doubly differential in the jet pr and rapidity. 

The charged particle multiplicity data from the analysis of minimum bias (MB) events ISJQ have 
already, and will continue to prove invaluable for the study of soft QCD and the underlying event 
(UE). Figure [T^ shows the azimuthal separation between the particles and the leading particle for MB 
events for different requirements on the leading particle pj. The increased coUimation of the event 
with increasing leading particle pr illustrating the onset of hard QCD can be clearly seen. For the 
Monte Carlo tune shown, the particle multiplicity transverse to the leading particle is too low and the 
MC events themselves appear more coUimated than the data. Figure[T}5 shows the inclusive single jet 
cross section doubly differential in the jet pj and rapidity. Within the large uncertainties, the NLO 
prediction corrected for hadronisation and UE describes the data well over five orders of magnitude. 

The LHC has been working well and after less than a year of collisions at 7 TeV is already 
providing a large range of valuable physics data. Despite the large range of high quality results 



On behalf of the ATLAS collaboration. 
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already available from the ATLAS Collaboration, these only begin to explore the available phase 
space. Many analyses are still statistically limited in the most interesting regions of phase space 
and the Collaboration is working hard to reduce the systematic uncertainties. Given the status of the 
statistical and systematic uncertainties, perturbative QCD appears to be in reasonable shape. 



First CMS results 



Thomas Hebbeke^ (RWTH, Aachen) 



In the year 2009 the LHC collider at CERN started with proton proton collisions at a center of 
mass energy of 900 GeV, later the energy was increased to 2.36 TeV and in March 2010 a 7 TeV 
run began which ended in November 2010. The CMS experiment Li3 | has recorded about 40/pb 
of integrated luminosity at this record energy. The CMS detector performed very well and many 
interesting measurements were made. Lead-lead collisions at a total center of mass energy of 574 TeV 
at the end of 2010 brought new insights into heavy ion physics. First of all the CMS collaboration 
has 're-discovered' all Standard Model particles, including W and Z bosons decaying leptonically, 
the heaviest quark top and the lighter quarks in form of various meson and baryon resonances. These 
analyses demonstrate that the detector has reached the design values for efficiency and resolution. The 
inclusive jet |T4J and dijet 1 15l production cross sections were among the first CMS measurements. 
The figure [14J shows that QCD calculations can reproduce the jet yields very well. 
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Figure 2: Jet spectra measured at various rapidities in pp at 7 TeV compared to NLO pQCD. 

The production of charged particles in minimum bias events (more precisely: in Non Single 
Diffractive events) was studied in great detail, in the pseudorapidity range |ri| < 2.4 and for transverse 
momenta as low as 30 MeV/c. In particular the yield as a function of r\, the multiplicity distribution 
and the pj distribution were determined at 0.9 TeV, 2.36 TeV and 7 TeV and compared to different 
model predictions and to data at other center of mass energies [161. Furthermore the production of 
strange particles like the were measured [ 17 |. Overall the increase of cross sections and multiplic- 
ity with center of mass energy is steeper than anticipated. The current models with parameters tuned 
without using LHC data do not provide a satisfactory description in all details. Also the shape of jets 



^On behalf of the CMS collaboration. 
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and the topology of hadronic events in general were analysed in detail II181I191 . The models PYTHIA 
and HERWIG provide a good description of these CMS measurements. A very interesting new feature 
was discovered in pp events with a very high charged particle multiplicity {N^^^ >110, ;?7' = 1 — 3 
GeV/c). In the two-particle correlation as a function of Ar] and Acj) a 'ridge' structure, a long range 
correlation in Ar] at small values of the azimuthal distance A(j) was revealed by CMS ll20l . Current 
Monte Carlo models cannot explain this feature. Finally already after a few days the first interesting 
heavy ion results were made public by CMS [21J: Z boson production in lead-lead collisions and 'jet 
quenching', seen as dijet events with very different energies carried by the two jets. 

Beyond QCD results, from the smooth falloff of the measured dijet cross section as a function of 
dijet mass one can set limits on new particles. For example excited quarks can be excluded within 
a contact interaction model up to a mass of 1.58 TeV at 95% confidence level |15|, thus improving 
older Tevatron limits significantly. From more than lOO'OOO W decays and about lO'OOO Z decays 
their production cross sections were measured, these results are in good agreement with NNLO QCD 
calculations [22]. Also a first measurement of the top cross section was made [23], confirming the 
expected strong rise (by a factor of about 25) with respect to proton - antiproton collisions at 2 TeV 
center of mass energy. 

Charged-particle multiplicity in pp and heavy-ion collisions at collider energies 

Jan Fiete Grosse-Oetringhaus ( CERN) 

The topical review ll24l summarizes and critically reviews measurements of charged-particle mul- 
tiplicity distributions and pseudorapidity densities in pp(^) collisions between ^/s = 23.6 GeV and 
1.8 TeV. Related theoretical concepts are briefly introduced: Feynman scaling which is based on phe- 
nomenological arguments about the exchange of quantum numbers predicts that the average number 
of charged particles increases with log^/s which implies that the rapidity density dNch/dy is a con- 
stant as function of ^/s. Feynman scaling is not fulfilled in the measured energy range. KNO scaling 
postulated by Koba, Nielsen and Olesen in 1972 asserts that the multiplicity distribution falls onto a 
universal curve when rewritten as P{Nch) —^1/ {Nch)P{z) with z = Nch/ {Nch)- KNO scaling is valid 
for NSD (non single diffractive) collisions in full phase space up to SppS energies {^/s = 200 GeV). 
For limited r] -intervals it still remains valid at LHC energies 111611251 . Negative binomial distributions 
(NBDs) describe multiplicity distributions well, a fact which is theoretically not well understood. 
Cluster models which assume an independent emission of clusters followed by stimulated particle 
emission within a cluster lead to multiplicity distributions of NED type. NBDs describe multiplicity 
distributions of NSD collisions in full phase space up to ^/s = 900 GeV where deviations have been 
observed. In limited rj-intervals the description still holds at LHC |'25','26 |. The combination of two 
NBDs, one representing soft and the other semi hard events (defined as events with and without mini- 
jets, respectively) describes multiplicity distributions up to the highest measured energy in full phase 
space (1.8 TeV) and in limited phase space intervals also at the LHC. The identification of trends of 
the fit parameters turns out to be ambiguous and assumptions are needed to obtain a coherent picture. 

Although no sound theory arguments exist at present why multiplicities in pp and e+e^ collisions 
should behave similarly, the average multiplicities as function of ^/s are similar when the concept 
of effective energy is introduced: E^s = ^/s — Eie^d.i — £^iead,2 where ^lead.i is the energy which is 
retained by the proton remnants. The available energy for particle production is characterized by the 
inelasticity K = Esf^/^/s. This concept allows one to fit the average multiplicity in pp collisions with 
the following form: fpp{\^) = fee{K^/s) + no. no is the contribution of the leading protons to the 
total multiplicity. One obtains a good fit result with K = 0.35 and no = 2.2, i.e., about one third of 
the energy is available for particle production compared to e"'"e~ collisions. Phenomenologically, one 
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could ask if this is a hint that only one out of the three valence quarks is available for particle produc- 
tion. However, the similarity between pp and e+e^ collisions cannot be found in more differential 
distributions, e.g. in (pseudo)rapidity densities. 

Results from LHC show a faster increase of the average multiplicity than anticipated by models 
112611271 . In particular phenomenological extrapolations and MCs with pre-LHC tunes underpredict 
the results at 7 TeV. Updated tunes which include LHC data indicate that parameters governing the 
amount of multiple-parton interactions need to be modified to reproduce the LHC multiplicities. 

LHCb QCD results 

Sebastian SchleicI^ (TU Dortmund) 

The Large Hadron Collider (LHC) delivered data in proton-proton collisions at unprecedented 
center of mass energies of ^/s = 900 GeV and 7 TeV, which allow one to test quantum chromodynam- 
ics (QCD) predictions in both, the perturbative and the non-perturbative regime. For example, the 
hadronization process falls into the latter. Predictions of the hadronization are based on phenomeno- 
logical models, that are mostly tuned on LEP data and their validity under LHC conditions needs to be 
confirmed by experiments. Designed for precision measurements in the B meson system, the Large 
Hadron Collider beauty (LHCb) experiment |28 | has several features that account for unique oppor- 
tunities for QCD studies in proton-proton collisions at the LHC: It covers a large forward rapidity 
range of 1.9 < r] < 4.9. Further, its tracking system includes a silicon tracker in close vicinity to the 
proton-proton interaction point and the experiment is equipped with a dedicated particle identification 
system based on two ring imaging Cerenkov detectors. 

The Ks production cross section measured at ^/s = 900 GeV in the kinematic region {px < 
1.6 GeV/c) x(2.5 <y < 4.0) is found |[29l to be in agreement with the Monte Carlo prediction 
(Pythia 6.4, Perugia tune ll30l . in the following referred to as MC), the pr spectrum tends to be 
slightly harder on data as compared to MC. Not as good agreement is found in inclusive (|) production 
cross section studied at ^/s = 1 TeV. The measured production cross section in the kinematical range 
(0.8 GeV/c <pt <5 GeV/c) x (2.44 <y< 4.06) is significantly enhanced with respect to MC, where 
also the pj spectrum is harder on data than on MC. 

The baryon suppression in A/ZTs is a sensitive test of fragmentation models because the initial 
state is purely baryonic. Similar considerations hold for particle-antiparticle ratios, since they probe 
the baryon transport from the beam to the final state. The particle ratio measurements K/Ks, A/A 
and p/p are presented at both, = 900 GeV and 7 TeV. The A/Ks ratio is underestimated by MC 
at both beam energies. In contrast, the A/A ratio in the kinematical range 2 < y < 4 is significantly 
overestimated by MC at ^/s = 900 GeV, whereas at 7 TeV the data is in better agreement, but still 
slightly lower than MC. The p/p ratio, measured in the range 2.8 < 3^ < 4.5, is slightly lower on 

= 900 GeV data than on MC, whereas it is in rather good agreement at 7 TeV. The hard QCD 
measurements presented are the Drell-Yan muon production pr spectrum at = 7 TeV, which is 
found to be in good agreement with the Monte Carlo prediction based on MCFM NLO. Additionally, 
the lepton px spectrum in in the W ^ fjv decay, as well as the charge asymmetry versus lepton 
pseudorapidity in this decay channel is presented. The results presented, still based on a relatively 
small data sample (59 nb^^), are in agreement with Monte Carlo predictions within their statistical 
uncertainties. 



On behalf of the LHCb collaboration. 
^AU results, except for the Ks cross section at ^/s = 900 GeV are preliminary 
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First TOTEM results and perspectives 



Karsten Egger^ (CERN & Cleveland) 



TOTEM is a dedicated experiment focused on forward physics complementary to the programmes 
of the large general-purpose experiments at the LHC OTTl . TOTEM will measure the total proton- 
proton cross-section with the luminosity independent method based on the Optical Theorem which re- 
quires a detailed study of the elastic scattering cross-section down to a squared four-momentum trans- 
fer \t\ of 10^^ GeV^ and the measurement of the total inelastic rate. Furthermore, TOTEM's physics 
programme aims at a deeper understanding of the proton structure by studying elastic scattering at 
large momentum transfers, and via a comprehensive menu of diffractive processes. To perform these 
measurements, TOTEM requires a good acceptance for particles produced at very small angles with 
respect to the beam. The coverage in the pseudorapidity-range of 3.1 < |ri| < 6.5 (rj = — lntan6/2) 
on both sides of the intersection point is accomplished by two telescopes for inelastically produced 
charged particles and complemented by Silicon detectors in special movable beam-pipe insertions - 
so called Roman pots (RP) - placed at about 147 m and 220 m from the interaction point, designed 
to detect leading protons at merely a few mm from the beam centre. 

During the year 2010, TOTEM has participated at the normal low-P* high-intensity runs with the 
vertical RP detectors at a distance to the beam of 18a and the horizontal at 20a. With this configura- 
tion, large-? elastic scattering could be measured for t- values above 2.2 GeV^. Furthermore, Single 
Diffraction and Double Pomeron processes were detected over a large region in the = Ap/p, t) 
plane and correlations between the forward proton and the particle densities in the very forward in- 
elastic detectors are studied as a function of ^. In addition, data were taken with only a few bunches 
during TOTEM dedicated runs where the vertical Roman Pots have been moved closer to the beams 
(7a) and a low intensity bunch with 10^" protons, to reduce pile-up, was added. Out of this data 
sample of several million events, 80k elastic scattering events with f-values above 0.4 GeV^ were 
extracted. The f-distribution showed the usual exponential slope at low t, but also exhibits a clear 
diffractive dip at around 0.6 GeV^, as it was observed for the first time at the ISR, whereas proton- 
antiproton elastic scattering only showed a shoulder in the f-distribution. The combination of these 
dedicated runs with the standard high-intensity runs will allow TOTEM to measure the f-distribution 
in the range of 0.4 - 4 GeV^ and to distinguish between the various models on the market. The low 
intensity bunch with a reduced pile-up of about 1% is used to measure the forward charged multiplic- 
ity and forward-backward multiplicity correlations. During the year 2011, TOTEM will concentrate 
on the total cross-section measurement, which becomes possible with the presently installed inelas- 
tic Cathode-Strip-Chambers. However, a special beam optics with large P* around 90 m has to be 
developed to enable the measurement of sufficiently small f-values necessary for the extrapolation 
to the optical point. Furthermore, the extensive studies of the forward particle flow and diffractive 
topologies will continue. 



LHCf is an LHC experiment designed to study the very forward production of neutral particles in 
pp collisions. Its results can provide valuable information for the calibration of the hadron interaction 
models used in Monte Carlo simulation codes, aiming in particular to clarify the interpretation of 

5 On behalf of the TOTEM collaboration. 
*On behalf of the LHCf collaboration. 



LHCf results 
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the energy spectrum and the composition of high energy cosmic rays as measured by air-shower 
experiments. The highest-energy data currently available for the forward neutral-pion production 
spectrum reach ,/s = 630 GeV (UA7 experiment [32] operated at the SppS in 1985-1986). The LHCf 
set-up consists of two imaging calorimeters (Arml and Arm2) symmetrically placed 140 m away on 
both the sides of the ATLAS interaction point, covering the pseudorapidity range |ri| > 8.4. Further 
information about the scientific goal, the technical details and the performance of the detectors can be 
found in the following references: |[33l439l . The experiment has successfully finished taking pp data 
at ^/s = 0.9 TeV and at ^/s = 7 TeV in 2009-2010. A preliminary analysis of the energy spectra of y- 
ray like and hadron like events measured by the calorimeters (Arml - upper plots, Arm2 - lower plots) 
at ^/s = 0.9 TeV, compared with the expectation of MC simulations, are shown in the Fig. [3l Only 
statistical errors are reported and simulations are normalized by total entries of y- and hadron-like 
events. The detectors were removed in July 2010 for an upgrade which will improve their radiation 
hardness, and they will be back in the LHC tunnel for the collisions at ^/s = 14 TeV. 
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Figure 3: Energy spectra of y- and hadron-like events measured by LHCf Arml (top) and Arm2 
(bottom) in pp at ^/s = 0.9 TeV compared with MC simulations. 



Tevatron results of relevance for cosmic-rays 

Lars Sonnenschei^ (RWTH Aachen) 

The two multi-purpose experiments D0 BOll and CDF f4ri are operated at the Tevatron collider, 
where proton anti -proton collisions take place at a centre of mass energy of 1.96 TeV in Run II. In 
the kinematic plane of 2^-scale and (anti-)proton momentum fraction x, Tevatron jet measurements 
cover a wide range, with phase space regions in common and beyond the HERA e7:)-collider reach. 
The kinematic limit of the Auger experiment is given by a centre of mass energy of about 100 TeV. 
Cosmic rays cover a large region of the kinematic phase space at low momenta x, corresponding to 
forward proton/diffractive physics and also at low scales, corresponding to the hadronisation scale and 
the underlying event. Therefore of particular interest are exclusive and diffractive measurements as 
well as underlying event, double parton scattering and minimum bias measurements. The kinematic 
limit of the Tevatron corresponds to the PeV energy region below the knee of the differential cosmic 
particle flux energy distribution. The data discussed here are in general corrected for detector effects, 

^On behalf of the CDF and D0 collaborations. 



10 



such as efficiency and acceptance. Therefore they can be used directly for testing and improving 
existing event generators and any future calculations/models. Comparisons take place at the hadronic 
final state (particle level). 

In particular for elastic and exclusive production measurements 1142 1143 1 forward proton detectors, 
which cover a pseudorapidity range of up to |ri| < 8 and momentum fractions of 0.03 < ^ < 0.10 
are useful, to detect the intact (anti-)proton. Many further analyses of exclusive and diffractive pro- 
duction ll44]450i have been accomplished. Underlying event, double parton scattering and minimum 
bias studies have been addressed by the measurements ||5TI - [56]| . The studies have been pioneering 
work in many cases. Methods have been established which are widely used by LHC experiments 
today. The measurements have provided very important input to theorists, in particular with respect 
to non-perturbative QCD physics, where phenomenological models are varying considerably. Most 
prominently the breakdown of factorisation between HERA and Tevatron has been established H6l . 
The double Pomeron exchange mechanism offers the possibility to study the exclusive Higgs produc- 
tion at the LHC, where predictions did vary by a factor of 1000 before the CDF measurement H6l . 
Already in Run I CDF has provided useful input for diffractive parton distribution functions. 

QCD results from HERA 

Armen Bunyatyai^ (Yerevan & DESY) 

New QCD results obtained by the HI and ZEUS experiments at HERA collider are reviewed. 
These results are based on data taken in e'^p collisions during 1994-2007 coiTcsponding to an inte- 
grated luminosity of almost 0.5 fb^^ for each experiment. HERA provides unique information on the 
proton structure. High center of mass energy ^/s = 320 GeV gives access to both the low Bjorken-jc 
domain and regime of high momentum transfers squared Q^. An ultimate precision of DIS cross 
section measurement is achieved by combining the HI and ZEUS measurements. The combined data 
are used as a sole input to a QCD fit to obtain HERAPDF sets. Important cross checks of the con- 
ventional QCD picture and the additional constraints for the gluon density distribution in the proton 
are provided by the measurements of the structure function and the charm and beauty cross sec- 
tions. New measurements of inclusive neutral and charged current scattering cross sections at high 
improve precision in this kinematic domain. Jet production at HERA provides an important testing 
ground for pQCD and new constraints for the gluon density distribution in the proton. The running 
of strong coupling is demonstrated and its value at mass, as{M'^), is determined from the jet 
measurements at HERA with high precision. 

The cross section of inclusive DIS diffractive process is measured within a wide kinematic range. 
The diffractive parton distribution functions (DPDFs) of the proton are determined from QCD fits to 
the data including the dijet production cross sections in diffractive DIS. Predictions based on these 
DPDFs are in agreement with the measured cross section of diffractive dijet and charm production in 
DIS at HERA and the longitudinal diffractive structure function F^^ . 

Data from the recent measurements of leading proton and neutron production are presented and 
compared to the theoretical models. The measurements are well described by the models which 
include the baryon production via virtual meson exchange. The leading proton and neutron data from 
the HI and ZEUS experiments are also compared with the hadron interaction models which are used 
in the analyses of ultra-high energy cosmic rays. The sensitivity of the HERA leading baryon data to 
the differences between the models is demonstrated. 



On behalf of the Hi and ZEUS collaborations. 
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Performance of the ALICE experiment for cosmic ray physics 

Bruno Alessandr^(INFN, Torino) 



A large number of atmospheric muon events were recorded during 2009 and 2010 for the cal- 
ibration, alignment and commissioning of most of the ALICE (A Large Ion ColUder Experiment 
at the CERN LHC) detectors. In this Workshop we presented the analysis of some of these data 
to understand the performances and the possibilities of ALICE to study topics connected to cosmic 
ray physics. The ALICE central detectors select atmospheric muons with zenith angle in the range 
0° — 60°. The muons are tracked in a large volume Time Projection Chamber (TPC) that measures 
the muon multiplicity, and for each muon the momentum, the sign, the direction and the spatial coor- 
dinates. An analysis of these observables and some correlation among them have been presented. In 
particular the muon multiplicity distribution and some events with very high multiplicity have been 
shown and detailed analyzed. A first attempt to measure the ratio / for vertical muons (0° — 20°) 
with a limited statistics has been presented and compared with world previous measurements. 

Horizontal muons, that is muons with zenith angle in the range 60° — 90° are very rare events that 
have been detected by the Forward Muon Spectrometer in 9 days of data taking. A selection of these 
events to obtain a good sample to measure the momentum distribution and the ratio /v+Z^^ at surface 
level has been discussed and measurements shown. 

CMS results on forward physics and other of relevance for cosmic rays 



Evidences of observation of single-diffraction at the LHC are presented at 900 and at 2360 GeV. 
Single-diffractive events appear as a peak at small values of the variable E±pz which is proportional 
to ^, the proton fractional energy loss, reflecting the 1 behaviour of the diffractive cross section. 
Single-diffractive events also appear as a peak in the energy distribution of the forward calorimeter 
HE, reflecting the presence of a rapidity gap over HE. The data have been compared on the detector 
level to PYTHIA 6 and PHOJET generators. PYTHIA 6 gives a better description of the non-diffractive 
component of the data, while PHOJET reproduces the diffractive contribution more accurately. 

The energy flow (at detector level) for minimum bias events and events having a hard scale de- 
fined by a dijet with Erjet >8 GeV (Erjet >20 GeV for ^/s = 1 TeV) in |ri| < 2.5 has been measured 
for the first time in hadron-hadron collisions in the forward region of 3.15 < |ri| < 4.9. The increase 
in energy flow in the forward region with increasing centre-of-mass energy is significant and is re- 
produced by Monte Carlo (MC) simulations for events with dijets, whereas it is underpredicted for 
minimum bias events. None of the MC simulations manages to describe all energy flow data. Monte 
Carlo tunes which are closer to the measurements for minimum bias events differ from measurement 
in dijet events. Particularly for minimum bias events, the measured energy flow a.t^/s = l TeV is larger 
than any of the predictions. Those MC simulations which best describe the energy flow in the forward 
region are different from those which best describe the complementary measurements of charged par- 
ticle spectra in the central region. The measurement of the energy flow in the forward region therefore 
provides further input to the tuning of MC event generators and it constrains modelling of multiple 
interactions at high energies. CMS measurements of charged particle pseudorapidity densities are be- 
ing compared with cosmic ray (CR) generators. The distribution proves of high discriminative power, 
revealing large spread of the CR MCs predictions at highest energy of 7 TeV. 

"^On behalf of the ALICE collaboration. 
"'On behalf of the CMS collaboration. 
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The flux ratio of positive- to negative-charge cosmic muons has been measured as a function of 
the muon momentum and its vertical component. This is most precise measurement below 100 GeV. 
The ratio was measured over broad range 10 GeV - 1 TeV of transition from approximately constant 
to a rising value. 

STAR results of relevance for cosmic rays 

Joanna KiryluI^(LBL, Berkeley) 

The Relativistic Heavy Ion Collider (RHIC) is a versatile accelerator situated at Brookhaven Na- 
tional Laboratory which commenced operations in 2000. It collides heavy ions as well as polarized 
protons at center of mass energies of up to 200 GeV per nucleon and 500 GeV, respectively. The 
Solenoid Tracker At RHIC (STAR) detector |[57l provides tracking, particle identification, and elec- 
tromagnetic calorimetry covering large acceptance. Key strengths include the capability to reconstruct 
jets, study correlations, and identify particles in high multiplicity environments. 

Measurements in unpolarized pp collisions at RHIC test perturbative QCD (pQCD) calculations. 
The data on jet and inclusive particle production cross sections II581I59II are in good agreements with 
Next-to-Leading-Order (NLO) pQCD calculations. Recent measurements of the boson produc- 
tion in 500 GeV pp collisions at RHIC and in 7 TeV collisions at the LHC are in good agreement with 
NLO pQCD over a wide range in [60]. 

The flux of prompt leptons at the Earth is of importance to cosmic ray (CR) and neutrino physics. 
Estimates depend strongly on models for the charm cross section and energy spectra. These models 
use the pQCD framework and extrapolate charm collider data to CR energies. The flux of prompt 
leptons is strongly dependent also on the small-jc nuclear gluon distributions. New dynamical effects, 
such as parton saturation that may be observed at forward collider rapidities, change the DGLAP dy- 
namics and thus the flux estimates 1611 . STAR and PHENIX data ll62l on electrons from heavy-flavor 
decays are consistent in the regions of kinematic overlap and are well described by Fixed-Order-Next- 
to-Leading Logarithm calculations |63|. The STAR J/*P data [64] are well -reproduced by calcula- 
tions ll65l using the color octet and singlet models in non-relativistic QCD. STAR has determined 
the B-hadron feed-down contribution to the inclusive J/^' yield from //^F-hadron azimuthal angle 
correlations. It is found to be 10-25% and has no significant ^/s dependence from RHIC to LHC 
energies [66|. The STAR T(lS-i-2S-i-3S) production cross section [,67,1 is consistent with world data 
and NLO pQCD calculations in the Color Evaporation Model ll68l . 

Collisions of dAu ions at ^/s^ = 200 GeV at RHIC have made it possible to study the modifica- 
tion of elementary QCD processes in cold nuclear matter (CNM), and provide insight in coherence 
effects or shadowing in nuclei, the saturation of small-x gluons, parton energy loss, and soft multiple 
scattering effects. Forward particle production is found to be suppressed [|59ll69l and back-to-back 
correlations are reduced llTOl . consistent with saturation models [ItTII . Hot matter effects have been 
studied at RHIC in AuAu collisions. The observation of phenomena such as jet quenching and collec- 
tive motion suggests that a thus far unobserved state of hot and dense matter with partonic degrees of 
freedom, resembling an ideal and strongly-coupled fluid, has been created [|72l . The ongoing Beam- 
Energy Scan program aims to observe the anticipated critical point in the QCD phase diagram ll73l . 
STAR has also observed anti-hypertriton production in AuAu collisions [|74l . the first ever observation 
of an anti-hypernucleus. The production and properties of antinuclei, and nuclei containing strange 
quarks, have implications spanning nuclear and particle physics, astrophysics, and cosmology. 



' 1 On behalf of the STAR collaboration. 
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2 Hadronic collisions at multi-TeV energies: Theory 



Monte Carlo tuning at the LHC 

Holser Schul¥^ (Inst, fur Phys., Humboldt-Univ. Berlin) 

Monte Carlo simulations of high energy physics processes are essential for many aspects of the 
LHC physics programme, e.g. the experiments use them to determine backgrounds to signal processes 
and to estimate reconstruction efficiencies, which are sources of systematic uncertainties that clearly 
dominate over statistical errors at the LHC. In order to reduce these systematic uncertainties, proton- 
proton collisions need to be simulated in such a way that they look as close to real data as possible. 
This however can in many cases only be achieved by optimising phenomenological model parameters 
to data, especially when 'soft' QCD effects such as multiple parton interactions or hadronisation are 
to be described. 

An overview of such parameter tuning strategies currently applied at the LHC has been presented, 
i.e. conventional manual parameter optimisation by the CMS collaboration and systematic tunings 
within ATLAS. The focus was clearly on the latter, where the software packages Professor and Rivet 
that allow for a systematic tuning effort using statistical techniques are used extensively. 

Rivet is an application that reads in generator independent events ("HepMC" format), processes 
these events by applying user written routines that mimic actual data analyses. After processing the 
generated events, histograms are produced that use the same binning as published data. 

Running a Monte Carlo event generator can be regarded as calculating a very expensive function. 
The key feature of Professor is the parameterisation of this expensive function by means of parameter- 
isations using polynomials such that one effectively produces a fast analytic model of the Monte Carlo 
generator response to shifts in parameter space. It is thus possible to get a very good approximate 
description of a generator at a certain point in parameter space in less than a second; a task that takes 
hours or even days with conventional methods. With this fast model, the task of tuning parameters 
is therefore passed on to constructing a goodness-of-fit measure between the parameterised generator 
response and real data. In this context, the necessity to have data, corrected for detector effects has 
been stressed. 

Further uses of the parameterisation like calculating the sensitivity of observables to shifts in pa- 
rameter space and an interactive Monte Carlo simulator ("prof-I") have been presented. The success- 
ful usage of Professor and Rivet has been illustrated by recent tunings performed within the ATLAS 
collaboration. Examples have been given ("AMBTl", "AUETl") and it has been stressed how fast 
the turn-around from taking new data to getting new tunings that include these data can be. Also, 
more special uses of Professor to study systematic variations ("Eigentunes", re-tuning using different 
PDFs) have been presented. 

Parton correlations and fluctuations 

Gosta Gustafso^^ ( Lund University) 

Multiple interactions and diffraction are important components in high energy collisions f 7511761. 
These effects are influenced by correlations and fluctuations in the parton evolution, and the under- 
standing of these features is therefore essential for a proper interpretation of data from LHC and 

'^On behalf of the Professor and Rivet collaborations. 

'■'in collaboration with C. Flensburg, L. Lonnblad, and A. Ster. 
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cosmic ray experiments. At high energies parton distribution functions at very small jc- values are 
governed by BFKL dynamics and saturation effects are important. Mueller's dipole model is a for- 
mulation of LL BFKL evolution in transverse coordinate space. The Lund Dipole Cascade model 
is a generalization of Mueller's model, which also includes non-leading effects from e.g. energy- 
momentum conservation and running coupling, saturation effects in the cascade evolution, and con- 
finement. The model is implemented in a MC called DiPSY, and in this talk I use it to study effects 
of correlations and fluctuations on double parton interactions and diffraction. 

In parton evolution a la BFKL the gluons are strongly correlated. An analysis of double parton 
distributions shows increased correlations for small x and large Q^. A spike develops for small sepa- 
rations between the partons in transverse coordinate space. The correlation can also be expressed in 
terms of an "effective cross section", which becomes reduced at high energies and large pj^. 

In the Good-Walker formalism diffractive excitation is determined by the fluctuations in the scat- 
tering amplitude between different components in the projectile wavefunction. In BFKL the proton 
substructure in terms of a parton cascade has large fluctuations and can fill a large rapidity range. 
An analysis of these fluctuations reproduces low and high mass diffractive excitation in DIS and pp 
collisions. For pp scattering the fluctuations are suppressed by unitarity constraints, which leads to a 
breaking of factorization between DIS and pp. 

The model can also be applied to nuclear collisions, and finally I present some preliminary results 
for exclusive final states in pp, pA, and AA collisions. 

A new model for minimum bias and the underlying event in sherpa 



Minimum bias events reveal not only the most complete view on the physics at hadron colliders, 
but also have an intimate connection to the underlying event and are thus highly relevant to many 
high-/7j^ processes. Higgs searches at the LHC, for instance, rely largely on event topologies with 
rapidity gaps. The feasibility of such measurements depends strongly on the survival probability of 
rapidity gaps. Apart from the connection to the underlying event, diffraction as an important part of 
minimum bias events is interesting in its own right. 

Unfortunately, no model describing soft, semi-hard, diffractive and hard QCD events has been 
implemented in a multi-purpose event generator so far. The Khoze-Martin-Ryskin model \J1\ is a 
multi-channel eikonal model that by summing all multi-pomeron diagrams is capable of describing 
elastic and inelastic scattering, low mass and high mass diffractive dissociation and central exclusive 
production. 

The Monte Carlo realisation relies on the partonic interpretation of the model. The simulation 
of elastic scattering is straight-forward, while the inelastic collisions are more involved. First, the 
number of exchanged ladders and the impact parameters of the ladders have to be generated. The 
emissions from the ladder are generated using a Sudakov form-factor, that accounts for absorptive 
corrections and Regge dynamics. The colour charge of the f-channel propagators has to be fixed, the 
singlet exchanges naturally give rise to rapidity gaps. Finally, the hardest emissions are corrected to 
pQCD matrix elements to reproduce the correct high energy behaviour. The model will be formulated 
also as a model for the underlying event and become available as part of the SHERPA 1.3 release. 



In collaboration with H. Hoeth, V. Khoze, F. Krauss, A. Martin and M. Ryskin 
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Multiple partonic interactions with herwig++ 



Stefan Gieseke (ITP, Karlsruhe) 



The focus of this talk is on the development of a multiple partonic interaction (MPI) model for 
minimum bias interactions and the underlying event in HERWIG++. We briefly summarize the general 
purpose Monte Carlo event generator HERWIG++ [ 78-81] before describing the development of the 
MPI model |[82l - [84l in detail. We explain the relevance of the main parameters of our model, the 
inverse radius and p™" in detail. The former characterizing the width of the spatial transverse 
parton distribution inside a hadronic projectile and the latter the transverse momentum down to which 
a partonic interaction will still be described by perturbative QCD. The development of the model was 
done in several steps which are briefly explained. 

1. A semi hard model for MPI, containing only interactions above /j™" II851I86II . 

2. A soft model, also for interactions below p™" ll87] - [89l . 

3. An extension of the model to allow for colour reconnections ||90l . 

The final model is shown to describe data from CDF ||9Tl[92l and recent non-diffractive min- 
imum bias and underlying event data from ATLAS at 900 GeV and 7TeV |l5l|7l[93l. The residual 
energy dependence of the model parameters is briefly discussed and an outlook to further work on 
this dependence is given. 

PYTHIA 8 status 

Torbjorn Sjostrand (Lund University) 

The Pythia 8 event generator iQ^I is the C++ successor to the Fortran-based PYTHIA 6 |[95l . 
frequently used in the study of pp/pp and e^e^ physics. One of the main developments is that 
Pythia 8 now contains a complete interleaving of multiparton interactions, initial-state radiation 
and final-state radiation, in one common sequence of decreasing p± scales. That is, the features at 
larger p± values set the stage for the subsequent dressing-up by softer emissions. The approximate 
matching of showers to hard-scattering matrix elements has been improved for a large set of hard 
processes (9^. For showers in QCD processes the first emission is compared with 2 — 3 matrix 
elements to confirm a reasonable rate 1971 . 

The traditional multiparton interactions framework is largely retained, but some new possibilities 
are added. One is that it is now possible to preselect two separate hard processes in the same event, to 
help simulate signals for double parton scattering. Another is that rescattering, where one parton scat- 
ters twice (or more) against partons from the other hadron, can now be simulated [98|. Unfortunately 
it is not simple to find a good experimental signal for such events. Other developments include an im- 
proved framework for the structure of diffractive events, a richer mix of underlying-event processes, 
and an updated set of parton distributions. 

Early attempts to tune Pythia 8 to minimum-bias data gave too much underlying-event activity. 
The problem has been traced to a double counting between some initial- and final-state radiation. This 
has now been fixed, and tunes to Tevatron data have been produced ll97ll . Unfortunately these tunes 
underestimate the activity observed in some of the early LHC data sets, which either may be owing 
to problems with the generator not reproducing different cut conditions, or point to some tension in 
the data. For now a slightly separate LHC tune has been made. 
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Pythia 8 does not reproduce the CMS ridge effect, which thus shows that some physics mech- 
anisms are still missing. Similarly there are problems e.g. with the particle composition and Bose- 
Einstein effects observed at LHC, that hints towards collective effects. 



Total hadronic cross sections at high energy 

Rohini Godbol^(CTS-Bangalore & CERN)) 

Energy dependence of total hadronic cross-sections is an important subject, both from a theoret- 
ical point of view due to its intimate relation to the non perturbative QCD dynamics and also from a 
phenomenological point of view, in the context of making accurate predictions for high energy cos- 
mic ray interactions based on the currently available information. In this talk I present a summary 
of the current state of data and model predictions for the same, paying particular attention to the 
Eikonalised Minijet Model (EMM) supplemented with soft gluon resummation ll99l . which tames the 
unacceptably strong energy rise of the EMM. 

The pPiPP and y/? (yy) cross sections, scaled by a VMD inspired factor of 330 ((330)^), all 
in fact show an almost universal behaviour, perhaps with a slightly faster rise for the yy induced 
processes flOO]. The models have to provide an explanation of the initial fall, normalisation at (and 
the position of) the minimum and the subsequent rise. The important issue to be addressed is the 
dynamics responsible for this rise, consistent with the Froissart bound. Then one can investigate 
the impact of these model predictions, extrapolated to cosmic ray energies. There exist different set 
of fits to the current data on total cross-sections, some of them with a form chosen so that various 
constraints from unitarity and analyticity are automatically satisfied. These have been then normally 
used to obtain the predicted total cross-sections at high energies. 

In the EMM models the rise of cross-section with energy is driven by the minijet cross-section 
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where 



calculated in perturbative QCD (pQCD), its rise with energy given by o,„inijet t— 

Pi III in 

£ = 7—1,7 being the degree of singularity of the gluon density in the proton. This has to be embedded 
in an eikonal formulation, which guarantees unitarity. The eikonalisation involves transverse parton 
overlap function in the two hadrons A&B, Aab{^) = f <i^Zjip^(Zji)pB(p — bi). This, along with the 
minijet (and some parametrisation of soft) cross-section, then is used in building the total cross- 
section, a^l^j^^.^ = 2/i/2g[i _f,-"(M/2]^ with n{b,s) = AAB{b,s)a{s). 

Different EMM models differ in the way AAB{b,s) is modeled. In the GGPS model f99), the 
energy dependent transverse space matter distribution is calculated as the Fourier Transform of the 
transverse momentum distribution of the partons, which is built through resummation of soft gluon 
emissions from the valence quark in the proton (to the leading order in as{Q^)) IllOli . In this BN 
EMM formulation, one of the important factor affecting the a'°' is behaviour of a, (2^) in the far 
infrared which is modeled by a form such that, a4(^^) — )• (1/^^)^ as — )• 0. The requirement 
that our form of be consistent with a confining potential and the singularity is integrable gives 
1/2 < p < \. Further, in the high energy limit, one can show a'°' ~ {elnsY^P and thus a high energy 
behaviour consistent with Froissart bound attains in our model naturally. Prediction of atot over 
the whole energy range requires also the soft cross-section, a*'°*'(5) and the corresponding overlap 
function A^g'*°^^ which are parameterised. 

The left panel of Fig. |4]taken from Ref. &9\ . shows comparisons of a variety of model predictions 
with data and each other, the blue band corresponding to the spread in our fits. We see that the LHC 



'^In collaboration with A. Grau, G. Pancheri and Y. Srivastava. 
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Figure 4: Predictions of the BN EMM (99) for pp/pp in the left panel and yp case 111 001 in the right 
panel, compared with different other model predictions as well as the data. 

data will already offer us a nontrivial discrimination among the models. For the photon induced 
processes, the BN EMM has to be supplemented by one more parameter Phad, the probability that 
the photon develops a hadronic structure. We find that the best fit seems to prefer a rise of yp cross- 
sections a little faster than the soft-pomeron predictions as well as those obtained, using factorisation, 
from our blue band for a'^^l of the left panel. Of course now the interesting thing is to try and see the 
effect of these range of predictions of o'^^J^'^ and the somewhat faster rise of the aj'^j on the cosmic 
ray simulations. Preliminary investigations II102I . for a photon energy of 10'^ eV, show that this 
difference can affect the development of the longitudinal profile of muonic showers by ~ 10%. 

The model has been recently applied to 71;? and nil cross-sections [103 j. In these cases the energy 
ranges over which the data are available are rather limited and in the latter case even the onset of 
the rise is not clearly established. The LHC with zero degree calorimeter (ZDC) can in fact provide 
nontrivial information to increase our confidence in making predictions for collisions with cosmic 
ray energies. In conclusion, the different models for total cross-sections seem to be good shape. 
The QCD based BN-EMM, which uses the experimentally measured parton densities and ideas of 
soft gluon resummation, in fact is able to even predict the Froissart bound. The same model gives a 
consistent description of high energy yp data and can be extended to n induced processes. The LHC 
data can play an important role in sharpening up our predictions so that they can be extended to the 
highest cosmic-ray energies. 

On strategies for determination and characterization of the underlying event 

Sebastian Sapeta (LPTHE, Paris) 

The underlying event (UE) is a soft activity which accompanies each hard process studied at 
hadron hadron collider. Its good understanding is of great importance since the UE affects a wide 
variety of high-z^r measurements, e.g. by introducing a bias or by degrading kinematic jet recon- 
struction. However, both an unambiguous definition of UE and its modeling faces a number of prob- 
lems. Therefore, it is particularly important to be able to measure the UE as well as possible so that 
corresponding experimental results can be used as an input to further constrain the models. 

We have carried out a twofold study devoted to this issue [104J- First, we asked the question of 
how the existing methods of UE determination perform on the practical aspects of the problem. De- 
veloping a simple toy model of UE and using is as a testing ground, we have examined two methods: 
the "traditional approach" [105] and the more recent area/median approach II106II1071 . One con- 
clusion from this part of our study is that for determinations of averaged quantities, like the average 
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transverse momentum of UE per unit area in the (y,(]))-plane, (p), both methods give comparably good 
results. In contrast, for event-by-event measurements and determinations of fluctuations of UE, the 
traditional approach is affected significantly more by the contamination coming from the hard part of 
the event. 

The second question we studied was that which observables related to the energy flow of UE 
are interesting to measure. Here, we chose the area/median method to examine more realistic UE 
from the Monte Carlo (MC) models and found noticeable differences between predictions of differ- 
ent generators/tunes extrapolated to LHC energy. Therefore, we conclude that a broader range of 
observables deserves dedicated measurements. Those include rapidity dependence of p, intra- and 
inter-event fluctuations and correlations. A first step of this program has already bin made by the 
CMS collaboration which used the area/median approach to measure the charged component of the 
UE at y/s = 0.9 TeV. The preliminary results 11108 II show that, even at low multiplicities, the method 
is capable to constraint MC tunes. 

Gluon saturation at the LHC from Color-Glass-Condensate 

Amir Rezaeian (Valparaiso) 

At high energy, a system of parton (gluons) forms a new state of matter: Color-Glass-Condensate 
(CGC) [ 109|. The CGC is the universal limit for the components of a hadron wavef unction which is 
highly coherent and extremely high-energy density ensemble of gluons. In the CGC picture, the den- 
sity of partons Pp with a typical transverse momenta less than Qs reaches a high value, «: I /a^ » 1 
(a^ is the strong coupling constant). The saturation scale is a new momentum scale that increases 
with energy. At high energies/small Bjorken-;c, Qs ^ where ^ is the scale of soft interaction. 
Therefore, {Qs) <C 1 and this fact allows us to treat this system on solid theoretical basis. On the 
other hand, even though the strong coupling becomes small due to the high density of partons, 
saturation effects, the fields interact strongly because of the classical coherence. This leads to a new 
regime of QCD with non-linear features which cannot be investigated in a more traditional pertur- 
bative approach. In the framework of the CGC approach the secondary hadrons are originated from 
the decay of gluon mini jets with the transverse momentum approximately equal to the saturation 
scale Qs{x) fllO]. The first stage of this process is rather under theoretical control and determines the 
main characteristics of the hadron production, especially as far as energy, rapidity and transverse mo- 
mentum dependence are concerned. The jet decay and hadronization unfortunately, could be treated 
mostly phenomenologically. 

The CGC [110] predicted 7 TeV data in pp collisions fT6l including (i) multiplicity distribution, 
(ii) inclusive charged-hadron transverse-momentum distribution, (iii) the position of peak in differ- 
ential yield 4) average transverse momentum of the produced hadron as a function of energy and 
hadron multiplicity. The same model also describes ep, eA and AA (at RHIC) data in an unified 
fashion supporting the universality of the saturation physics. It has been shown that the observed 
ridge phenomenon in pp collisions at the LHC can be also explained by the CGC [111]. There exists 
some ideas how to simulate the CGC state in AA collisions due to higher density but evidence for the 
formation of the CGC state (gluon saturation) in proton-proton interaction will be a triumph of the 
high-density QCD and the CGC. 

The physics of AA collisions is more complicated compared to pp{A) and ep{A) collisions. The 
ALICE collaboration has recently released new data for the multiplicity in central PbPb collisions at 
y/s^ = 2.76 TeV 111121 . There are some surprises in the ALICE data: (i) the power-law behavior 
on energy in AA is so different from pp collisions which is not very easy to accommodate within 
the CGC approach, (ii) the models that describe DIS for proton, DIS for nucleus, the LHC data 
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for proton and RHIC data apparently failed to describe the ALICE data with the same accuracy. It 
appears apparently to be difficult to describe at the same time, HERA and RHIC and the new ALICE 
data for the multiplicity [ 1 13 1. First notice that, the ALICE — 5% centrality bin at ^JJ^ = 2.76 TeV 
corresponds to Npart = 381 while our approach based on the Glauber model gives Npart = 374 HI 141 . 
Therefore, our actual prediction llll4t for the same centrality bin will be higher than what the ALICE 
collaboration quoted in their paper II112I . Assuming that the ALICE data is correct, saturation models 
gave correct predictions for multiplicity in AA collisions at the LHC within about less than 20% 
error. Indeed, this is not horribly bad given the simplicity of the approach. However, this will give 
rise to several open questions: (i) what is the role of final-state effects? (ii) how the mini-jet mas 
changes with energy/rapidity in a very dense medium?, (iii) what is the effects of fluctuations and 
pre-hadronization?. One should also have in mind that the kj factorization for AA collisions has 
not yet been proven and gluon production in AA collisions is still an open problem in the CGC. To 
conclude, the first LHC data for AA collisions has already created much excitement in the heavy ion 
community and it opened a fresh and hot debate on how the saturation physics changes from ep, pp 
and eA to AA collisions. 

Systematic study of inclusive hadron production spectra in collider experiments 

Audrey Rostovtse^ (ITEP, Moscow) 

There exists a large body of experimental data on hadron production in high energy proton- 
(anti)proton, photon-proton, photon-photon and heavy ion collisions. In the present report the ex- 
perimentally measured inclusive spectra of long-lived charged particles produced at central rapidities 
in the colliding particles center of mass system are considered. The analysed published data have 
been taken with a minimum bias trigger conditions and at center of mass energy (y^) ranging from 
23 to 2360 GeV. 

The charged particle spectra as function of transverse momentum are traditionally approximated 
using the Tsallis-type (power law) function. However, a closer look at the fits to the available data 
discloses systematic defects in this approximation. It is found, the parameterization 

=A,exp(-4-/r,) + (1) 

PrCipT + 

is in much better agreement with the data then the Tsallis approximation. The variables in the equation 
above are self-explanatory. The most surprising feature of the new parameterization ([Til is a strong 
coiTclation between the parameters Tg and T. Though the physical origin of the observed correlation is 
not quite clear, it provides an additional constraint for the parameterization ([T]) and therefore reduces a 
number of free parameters. Interestingly, a similar combination of the Boltzmann-like and power-law 
terms is observed in the photon energy spectra from the sun flares. 

The relative contributions of the terms in ([Hi are characterized by a ratio R of the exponential 
to power law terms integrated over pj-. Interestingly, for pp and pp data this ratio R is almost in- 
dependent of the collision energy and equals to about 4, while for AuAu it reaches minimum values 
(about 2) at medium centralities of heavy ion collisions. In addition, in the high energy DIS, photo- 
production and yy collisions the power law term of the new proposed parameterization ([Tjl dominates 
the produced particle spectra. Thus, only the inclusive spectra of charged particles produced in pure 
baryonic collisions require a substantial contribution of the Boltzmann-hke exponential term. 

'^In collaboration with A. Bylinkin. 



20 



Finally, a map of the parameters T and n for proton-(anti)proton, heavy ion, y-proton and yy 
collision at different energies is drawn. There are two clearly distinct trends seen on the map. The 
pp and pp collision data show an increase of the T— parameter and decrease of the n parameter with 
collision energy increasing. The second trend, where the values of both parameters the T and n 
increase, is defined mainly by the RHIC AuAu collision data at -.Js^ = 200 GeV per nucleon. In this 
case a simultaneous increase of the T and n values corresponds to an increase of the centrality of heavy 
ion collisions. Surprisingly, the both trends cross each other at medium centralities corresponding to 
the minimum bias AuAu collisions and p'p interactions with energy of y/S^=200 GeV. Naively one 
could expect the single p'p interaction has more similarity to the very peripheral single nucleon- 
nucleon interactions. Contrary to that, DIS, y;? and yy interactions belong to the second trend and 
are located on the parameter map nearby very peripheral heavy ion interactions at about the same 
collision energy per nucleon. A more extended version of this report can be found in II115II . 

Hyperon transverse momentum distributions m pp and pp collisions 

Olga Piskounova (Lebedev Inst., Moscow) 

The analysis of data on hyperon transverse momentum distributions, dN jdpj, that were gathered 
from various experiments (WA89, ISR, STAR, UAl and CDF) reveals an important difference in the 
dynamics of multiparticle production in proton-proton vs. antiproton-proton collisions in the region 
0.3 GeV/c < pt <3 GeV/c. Hyperons produced with proton beams display a sharp exponential slope 
at low Pt, while those produced with antiproton beam do not. Since LHC experiments have proton 
projectiles, the spectra of multiparticle production at the LHC 1271 should be "softer" in comparison 
to predictions, because the MC predictions were based on Tevatron (antiproton) data. 
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Figure 5: Strange baryon transverse momentum distributions from different experiments. 

The available data of many high energy experiments on pp colUsions Illl6[|117 l as well as on pp 
collisions of lower energies 1 118111191 and neutron-carbon reaction ri201 are considered in this article 
in order to understand the influence of quark composition of beam particle on the shape of transverse 
momentum spectra of A" hyperon production as at high energy collider experiments as in low energy 
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fixed-target experiment (Fig. The difference in pj- spectra of A'^'s produced in higli energy pp and 
pp collisions can not be explained in the QCD theoretical models, because at the collider energies 
both interactions should give the multiple particle production due to Pomeron (or multi-Pomeron) 
exchange. The total cross section and the spectra in pp and pp collisions are to be similar because 
of the Pomeron exchange between two interacting hadrons that should not be sensitive to the quark 
contents of colliding beams at high energies. 

Unfortunately, this difference was not studied enough at ISR, where both projectiles were avail- 
able. The important fact is that the latest collider experiments were carried out with antiproton beams. 
It was mistake to suggest that pp and pp at high energy are giving the similar transverse momentum 
distributions. The spectra of hyperons that are produced with proton beam have a sharp exponential 
slope at low pj, while the spectra with antiproton beam have not. The highest energy experiments 
(UAl and Tevatron) shows harder pj spectra, that is not only the result of growing energy - it is the 
result of different form of transverse momentum distributions in different reactions. 

From the point of view of the Quark-Gluon String Model (QGSM) [il2H . the most important con- 
tribution to particle production spectra in antiproton-proton reactions is due to antidiquark-diquark 
string fragmentation. Baryon hadroproduction spectra are sensitive to quark-diquark structure of in- 
teracting hadrons as well as to the energy splitting between these components. Asymmetric reactions 
may provide us with a new "stereoscopic" view on the hadroproduction mechanism. Measurements of 
Pt spectra in antiproton-proton interactions at a variety of energies can thus constrain the contribution 
from the fragmentation of antidiquark-diquark string. This study may have impact not only on the in- 
terpretation of LHC results, but also on cosmic ray physics and astrophysics, where matter-antimatter 
asymmetry is being studied. 

DPMJET-III and pp data from the LHC 

Johannes Ranft (Siegen University) 

Monte Carlo codes based on the two-component Dual Parton Model (soft hadronic chains and 
hard hadronic collisions) are available since 10-15 years: PHOJET for hadron-hadron Qi-h) and 
photon-hadron (y-/?) coUisions II122II and DPMJET-III based on PHOJET for hadron-nucleus Qi-A) 
and nucleus-nucleus (A-A) collisions II123II . 

At the LHC particle production m pp collisions was measured by the three Collaborations CMS, 
ALICE and ATLAS. Here we compare with measured pseudorapidity distributions at 900, 2360 and 
7000 GeV cm. energy, with pj distributions of charged hadrons at the same energies and with p 
to p ratios. The problem to be solved at the beginning was, that DPMJET-III and all other event 
generators did predict pseudorapidity distributions rising slower with energy than the LHC data. To 
solve this problem we had to redetermine the parameters of DPMJET-III in such a way, that agreement 
with the data is achieved. We can present at the moment only one preliminary solution. In this 
solution we introduce an energy dependence in two of the DPMJET-III parameters. We call this a 
preliminary solution, since we think that also the new parameters should not depend on the energy. 
We are confident, that we will find soon such a solution, but at the moment all solutions with energy 
independent parameters have still problems. 

We were able to present at the meeting charged pseudorapidity distributions which agree perfectly 
with the CMS data. We found also central antiproton to proton ratios in agreement with the ALICE 
data. The pj distributions of DPMJET-III agree with the CMS distributions. But one weakness of the 
Pt comparisons is, that at present the average pt values at energies lower than the LHC energies are 
slightly higher than the data, this is also one problem which has to be solved. Finally, we find perfect 
agreement of charged hadron multiplicity distributions comparing DPMJET-III with the ALICE data. 
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3 Cosmic-rays at Ultra-High Energies: Experiments 



Results from the Pierre Auger Observatory 

Lorenzo Caz.oi^ (LIP. Lisbon) 

The Pierre Auger Observatory is the largest cosmic ray (CR) observatory on Earth, covering 3000 
km^ of the high plateau in the Argentinian region of Pampa Amarilla II124L It detects CR air show- 
ers in two complementary ways: an array of water-Cerenkov tanks samples the secondary particles 
at ground and fluorescence telescopes observe the longitudinal development of the electromagnetic 
cascade |[T25l . 

An energy spectrum has been recently published HI 261 covering the energy range from 10^^ eV 
to above lO^'' eV. The dominant systematic uncertainty stems from the overall energy scale, and 
is estimated to be 22%. The position of the ankle at logiQ{Eankie/^V) = 18.61 ±0.01 and a flux 
suppression above log (£1/2/^^) = 19.61 ±0.03 have been determined. The suppression is similar 
to what is expected from the GZK effect for protons or nuclei as heavy as iron, but could also be 
related to a change of the injection spectrum at the sources. 
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Figure 6: (Xmax) and RMS{X,„ax) measured by Auger compared to different hadronic interaction 
models. 

Analyses of both the mean values and the fluctuations of the shower maximum (X^^v) (plotted 
in Fig.|6l in comparison with the prediction of different hadronic interaction models) reveal a change 
in the energy dependence of the composition around the ankle and a gradual increase of the average 
mass of cosmic rays with energy, provided that there are not significant changes in the properties of 
the hadronic interactions at ultrahigh energies. More details can be found in P 127II . 

A detailed comparison of the muonic and electromagnetic content of the air showers far from the 
core shows an excess in the number of muons of ~ 50% compared to protons simulated with QGSJET- 
II (~ 30% excess compared to iron) III281 . This analysis also suggests a larger energy scale (27%) 
with respect to the fluorescence detector, being still compatible within the systematic uncertainties. 

Showers initiated by photons and neutrinos have distinct signatures compared to showers initiated 
by protons or other nuclei, being possible to discriminate them ||129[ll30l and set bounds on their 
fluxes. The neutrino and photon bounds already exclude 'top-down' models for the production of 
ultra-high energy cosmic rays, favoring acceleration in astrophysical scenarios. 

In 2007, it was shown that the arrival directions of CRs with energy in excess of 55 EeV corre- 
lated with the positions of nearby AGN 111311 . As further data has been added, the degree of corre- 

l^On behalf of the PAO Collaboration. Support FCT-Portugal (CERN/FP/109286/2009) and ECT* are acknowledged. 
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lation has decreased from 69_i^% to 38_^%, to be compared with the 21% expected if the flux were 
isotropic 1 132 1. The region of the sky with the largest observed excess with respect to isotropy corre- 
sponds to Cen A, which is the closest AGN. Nevertheless, at present there ai^e multiple astrophysical 
models of anisotropy which are fully consistent with the observed distribution of arrival directions. 
While a correlation of arrival directions with nearby matter on small angular scales is plausible for 
protons above 55 EeV, it is puzzling if the CRs are heavy nuclei as suggested by the X,„ax mea- 
surements, since they are expected to undergo large deflections due to the galactic and extragalactic 
magnetic fields. Definitive conclusions must await additional data. 

HiRes results: The final word (almost) 

Douglas R. Bersmari^ (Univ. of Utah) 

The High Resolution Fly's Eye (HiRes) experiment was a ultra-high energy cosmic ray detector 
operating in Utah from 1997-2006. It used the fluorescence technique to detect cosmic ray air showers 
in stereo using two detector sites. 

HiRes was the first experiment to observe the GZK cutoff. We observed the break with a signif- 
icance of greater than 5a in monocular mode 13 and confirmed that discovery in stereo mode 11331 . 
The position of the break is consistent with what is expected for extragalactic protons as measured 
with Beresinki's £1/2 test 111341 . 

HiRes has made a direct measurement of cosmic ray composition at energies above 10^^^ eV 
using both the average Xmax technique and by measuring the width of the X^ax distribution 11351 . 
Both show a composition consistent with pure protons. We note that great care must be taken to 
account for both trigger and reconstruction biases for the mean X^ax measurement. We also note 
that the RMS is a biased estimator of the width of a non-Gaussian distribution; we use a Gaussian fit 
to the central part of the distribution (within twice the RMS of the mean) as a better estimator of the 
width. HiRes also measured the average width of individual showers. This also agrees with a protonic 
composition. 

HiRes has made a number of searches for anisotropy in the arrival directions of cosmic rays. We 
observed none: no correlation with AGN's II136I . no correlation with the local large scale structure as 
determined by galaxy surveys II137I . 

Telescope-Array results 

Nobuyuki Sakura^ (Osaka City Univ.) 

Telescope Array (TA) is a largest ultra high energy cosmic ray detector in northern hemisphere. 
Operation was started about 3 years ago, and the gathered data have delivered the interesting results 
on the extremely high energy cosmic rays. The detector consists of 507 plastic scintillation detectors 
(SD) which cover the ground area of 680 km^ in 1.2 km mesh and 3 fluorescence telescope stations 
(FD) which surround the scintillator detector array and look inward. Each detector of SD has 2 layers 
of plastic scintillator plate of 3 m^ area and 1.2 cm thickness. For FD, Telescope Array adopts two 
different types of telescopes. Two FD stations, which are called as BRM-FD and LR-FD, are newly 
developed for Telescope Array. One FD station which is called as MD-FD consists of the telescopes 
which had been operated as Hires-I detector. 



'^On behalf of the HiRes collaboration. 

^On behalf of the Telescope-Array collaboration. 
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We analyze energy spectra of ultra high energy cosmic ray using 3 different data sets: FD-mono, 
Hybrid and SD. In the FD-mono analysis, MD-FD data is analyzed by the same program as HiRes-I. 
And we obtained energy spectrum which is consistent with HiRes result. This means that HiRes-I 
detectors which was moved to Telescope Array site have reproduced HiRes result. The Hybrid energy 
spectrum is obtained from the analysis which use both of FD data and SD data in order to improve the 
geometrical reconstruction, although its energy is reconstructed as same as FD analysis. The energy 
spectrum obtained by the hybrid analysis is consistent with FD-mono result. TA SD energy spectrum 
is also consistent with spectra from both of FD-mono analysis and the hybrid analysis. In analysis 
of SD data, SD energy scale is scaled so as to agree with FD energy scale using the result of hybrid 
analysis. Primary composition is studied using the shower maximum {X,„ax) observed by FD stereo 
data. In the energy region of 10^^-^ ~ 10'^-^ eV, the averaged X,„ax is consistent with the proton 
primary hypothesis. Arrival direction study shows no correlation with Active Galactic Nuclei so far. 
The existence of ultra high energy photon is studied using the shower front curvature observed by SD, 
but no candidate is found in data. Electron LINAC (ELS) is installed in front of BRM-FD to calibrate 
using the electron beam in this year. The ELS calibration is expected to improve the systematic error 
of FD drastically. 

ARGO-YBJ results 

Ivan De Mitr^ (INFN, Lecce) 

Cosmic ray physics in the 10^^ — 10^^ eV primary energy range is among the main scientific goals 
of the ARGO-YBJ experiment f l381ll39l . The detector, located at the Cosmic Ray Observatory of 
Yangbajing (Tibet, RR. of China) at 4300 m a.s.l., is a full coverage Extensive Air Shower array con- 
sisting of a carpet of Resistive Plate Chambers of about 6000 m^. The apparatus layout, performance 
and location offer a unique possibility to make a deep study of several characteristics of the hadronic 
component of the cosmic ray flux in an energy window marked by the transition from direct to indi- 
rect measurements. In this short summary we will focus on hadronic interaction studies that are being 
performed within the experiment in the primary energy range going from 1 TeV to 1 PeV. 

The proton-air cross section has been measured. The total proton-proton cross section has then 
been estimated at center of mass energies between 70 and 500 GeV, where no accelerator data are 
currently available. Other hadronic interaction studies can be performed by exploiting the detector 
capability to have very detailed information on the shower front space-time structure and the lateral 
distribution function by also using the analog readout of the RPC's. Because of lack of space, here 
we will report on the p-air and pp cross section measurement only. 

The measurement is based on the shower flux attenuation for different zenith angles, i.e. atmo- 
spheric depths 1 140 1 . The detector location (i.e. small atmospheric depth) and features (full coverage, 
angular resolution, fine granularity, etc.) ensure the capability of reconstructing showers in a very 
detailed way. These features have been used to fix the energy ranges and to constrain the shower 
ages. In particular, different hit (i.e. strip) multiplicity intervals have been used to select showers 
corresponding to different primary energies. At the same time the information on particle density, 
lateral profile and shower front extension have been used to select showers having their maximum 
development within a given distance/grammage Xj„, from the detection level. This made possible the 
unbiased observation of the expected exponential falling of shower intensities as a function of the 
atmospheric depth through the sec 6 distribution. After the event selection, the fit to this distribution 
with an exponential law gives the slope value a, connected to the characteristic length A through the 

2°On behalf of the ARGO-YBJ collaboration. 
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relation a = /tq/A. That is: 

/(e) = A(e) /(e = o) e-^o^^^^-i) (2) 

where A(0) accounts for the geometrical acceptance of each angular bin. The parameter A is con- 
nected to the proton interaction length by the relation A = kXint, where k depends on hadronic inter- 
actions and on the shower development in the atmosphere and its fluctuations II141II . The actual value 
of k must be evaluated by a full MC simulation and it depends also on the experimental approach, 
the primary energy range and on the detector response. The p-div production [142] cross section 
is then obtained from the relation: ap_^,v(mb) = 2.41 x 10^/?i/„^(g/cm^), while several theoretical 
approaches can be used to get the corresponding pp total cross section Opp II143II . 
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Figure 7: The pp total cross section obtained by ARGO-YBJ, together with results published by other 
cosmic ray and accelerator experiments II140II . 

Results are shown in Fig. |7] As can be seen, ARGO-YBJ data lie in an energy region not yet 
reached by pp colliders (and still unexplored by p^ experiments), favouring the asymptotic ln?-{s) 



increase of total hadronic cross sections as obtained in II144II from a global analysis of accelerator 
data. 



KASCADE-Grande results 



Paul Dol^ (KIT, Karlsruhe) 



Testing of hadronic interaction models QGSJET-II-2 and EPOS 1.99 implemented in the CORSIKA 
program have been performed with KASCADE-Grande air shower data in the energy range of lO'^ 
to 10^^ eV II1451I1461 . From the muon density investigations, the EPOS 1.99 model indicates that 
light abundances of primary cosmic ray particles would be needed to fit the data. On the other hand, 
the QGSJET-II-2 model describes the data with an intermediate primary abundance between proton 
and iron nuclei. The reconstructed all-particle energy spectra are presented by using the hadronic 
interaction models QGSJET-II-2 and EPOS 1.99. The resulting spectra show that the interpretation 
of the KASCADE-Grande data with EPOS 1.99 leads to significantly higher flux as compared to the 
QGSJET-II-2 result. More detailed investigations of EPOS 1.99 are still in work. 



On behalf of the KASCADE-Grande collaboration. 
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IceCube results 



Lisa Gerhard^ (LBL, Berkeley) 



High energy neutrinos offer a unique view of distant, energetic astrophysical objects, as they are 
neither bent by ambient magnetic fields nor absorbed by the interstellar medium. Possible sources 
of neutrinos include active galactic nuclei, gamma ray bursts, the highest energy cosmic rays, and 
interactions of exotic objects. The IceCube neutrino detector uses the ice at the South Pole as a 
Cerenkov medium for the detection of high energy neutrinos. It is composed of an in-ice, three- 
dimensional array of photomultiplier tubes [147] and a surface air shower array. Construction of 
the IceCube detector began in 2005 and was finished in 2010, bringing the detector to its full cubic 
kilometer size. Using data from the partially constructed detector, the IceCube Collaboration has 
searched for point sources of neutrinos II148II149I and found results consistent with the expectation 
from the background of atmospheric neutrinos. It has set stringent upper limits on the diffuse fluxes 
of extremely high energy neutrinos [150] and on the flux of neutrinos in coincidence with Gamma- 
Ray Bursts II151II and has set limits on the accumulation of dark matter in the Sun II1521I1531 . It has 
measured the flux of atmospheric neutrinos up to 400 TeV 111541 and the anisotropy of the arrival 
directions of cosmic rays with a median energy of 20 TeV II155II . IceCube data collection continues. 



On behalf of the IceCube collaboration. 
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4 Cosmic-rays at Ultra-High Energies: Theory 



Open problems in cosmic ray pliysics and tlie importance of understanding liadronic interac- 
tions 

Ralph Engel (KIT, Karlsruhe) 

In Fig. m a compilation of measurements of the all-particle spectrum of cosmic rays is shown 
(from fl561 . updated). The most striking features are the knee at about 3 x 10^^ eV, the ankle be- 
tween 10^^ - lO^^eV, and the suppression of the flux at the very highest energies. Understanding the 
origin of these characteristic breaks in the power-law of the flux is key to identifying the galactic 
and extragalactic sources of cosmic rays and the corresponding particle acceleration and propagation 
mechanisms. 
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Figure 8: The all-particle flux of cosmic rays. For references to the data see II156II . 

For example, if the knee in the flux stems from features of diffusive shock acceleration processes 
it is expected that the fluxes of individual elements will exhibit knee-like features shifted in energy 
according to the magnetic rigidity ~ E'/Z. Similarly, particle physics scenarios lead to the prediction 
of a scaling with mass number ~ E/A. Many other models have been developed for the knee, for 
example 1157111581 - for a review, see fT59]. Similarly, there are several interpretations of the ankle. 
It seems natural to assign this feature to the transition between galactic and extragalactic cosmic ray 
sources II160II . Alternatively, in the dip model the ankle is the imprint of e+e pair production in 
extragalactic propagation II161I . requiring a proton-dominated composition. 

All these model scenarios differ in the predicted composition, making the accurate measurement 
of the mass composition and its evolution with energy a prerequisite to making progress in the field. 
Currently the largest uncertainty in the composition interpretation of air shower data is related to 
the description of hadronic interactions that has to be done with phenomenological models II162I . 
Measurement of multiparticle production at HERA, RHIC and LHC and using the data for refining 
existing interaction models will allow us to make significant progress in reducing this uncertainty. 
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UHECRs and hadronic interactions 



Paolo Lipari (Roma) 



The interpretation of the data on the Uhra High Energy Cosmic Rays requires an understanding of 
the development of hadronic showers and therefore a sufficiently accurate description of the properties 
of hadronic interactions. Since the energy spectrum of CR extends up to £" ~ 10^" eV, this requires 
an extrapolation to the cm. energy oi ^/s ^ 430 TeV. 

The observations of fluorescence light, pioneered by the Fly's Eye collaboration, allows one to 
estimate the energy of a CR particle in a quasi-calorimetric way, with only little model dependence. 
The study of the shape of the longitudinal shower development allows one in principle to determine 
the mass number A of the primary particle, but is also strongly dependent from the hadronic interaction 
properties. The position of the maximum of the shower longitudinal development ^max 

is a good 

indicator of A. The average Xmax for a primary particle of energy E and mass A, in reasonably good 
approximation, is given by: 



(the quantity Dp is known as the "elongation rate"). This equation allows one to estimate the average 
mass (or (logA)) and the evolution with energy of the composition of UHECR, that is J(logA) /dE, if 
one has a good theoretical control of the model dependent quantities Xq and Dp. It is therefore neces- 
sary to estimate the systematic uncertainty in the calculation of these quantities. It is also interesting 
to discuss what properties of the hadronic interactions determine Xq and Dp. The main contributions 
come from the hadron interactions lengths (including the interaction lengths of mesons) and the inclu- 
sive energy spectrum of secondary particle in the projectile fragmentation region. The results that can 
be obtained at LHC with a precision measurement of the pp cross section and of particle production 
properties will be important in constraining the models. 

It is interesting to note that it is in principle possible to determine the composition of CR without 
any use of shower development, and therefore use the CR measurements to obtain information about 
the properties of hadronic interactions. One possibility is the observation of the imprints of energy 
losses on the observed energy spectrum (since the kinematical thresholds for the processes of e^e~ 
pairs and pion production are A dependent), and the estimate of the deviation due to astronomical 
magnetic fields (that depends on Z). This could in principle allow one to determine properties of 
hadronic interactions from CR observations. The observations are at this point inconclusive. The 
energy scale of the HiRes detector is consistent with the hypothesis of attributing the "ankle" spectral 
feature, and the high energy suppression with energy loss imprints on a smooth spectrum strongly 
dominated by protons, but this interpretation is not consistent with the energy scale of the AUGER 
experiment. The AUGER collaboration has observed a correlation between the direction of the high- 
est energy particles with potential extragalactic sources but the interpretation of the results remains 
ambiguous. 

The study of the width of the distribution in the measurement of Zmax also allows one to estimate 
the mass composition of CR, since the development of large A primary particle has smaller fluctu- 
ations. The results of the AUGER experiment suggest that the highest energy particles are large A 
nuclei, however these results are not confirmed by results of the HiRes and Telescope Array collabo- 
rations. 




-Zo+Dp log - =XQ + Dp \ogE-Dp\ogA 




(3) 
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On the relation between air shower predictions and features of hadronic interactions 

RalfUlrich (KIT, Karlsruhe) 



The nature of cosmic ray particles at the very highest energies is still not understood. Even with 
experiments like the Pierre Auger Observatory fl63], HiRes [164] and Telescope Array fl651 de- 
livering large quantities of high quality data, it is not straightforward to interpret these data. The 
fundamental problem is that the observations are very indirect: Only extensive air shower cascades, 
which are initiated by the cosmic ray primaries, are observed. For an accurate analysis of these air 
shower data a detailed understanding of interactions in the cascades is required. The particle produc- 
tion characteristics that are important in the context of air showers are interactions at energies of up 
to ^/s ~ 350 TeV and particle production in the very forward direction (rj > 6). Particle production 
characteristics in this phase space have a strong impact on the modelling of the evolution of air shower 
cascades 111661 . 
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Figure 9: Measured average depth of the shower maximum 11167111681 of UHE cosmic rays compared 
to SIBYLL predictions with a modified extrapolation of the elasticity to ultra-high energies. 

In Figure |9] we demonstrate this at the example of the elasticity, ^eia = E^ax/Etou of interactions. 
It is shown that the interpretation of cosmic ray data relies on the detailed understanding of hadronic 
interaction physics. The LHC has the potential to significantly reduce the uncertainties because of 
two reasons: Firstly, it operates at energies that are already very relevant in terms of cosmic ray obser- 
vations and, secondly, it has significant capabilities to study the forward phase space of multiparticle 
production that are relevant for the air shower modelling. 

The event generator sibyll 2.1 

Ralph (KIT, Karlsruhe) 

SIBYLL is an event generator optimized for simulating high energy interactions needed for the 
description of extensive air showers and for the calculation of inclusive muon and neutrino fluxes. In 
comparison to other event generators, a rather basic and straightforward model for the description of 
multiparticle production is implemented. The initial version of the model 1 1691 was upgraded in 2000 

■^^In collaboration with Eun-Joo Ahn, Thomas K. Gaisser, Paolo Lipari, and Todor Stanev. 
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to include post-HERA parton density parametrizations and to introduce multiple soft interactions and 
a better treatment of diffraction dissociation 1 170|. 

Proton/pion/kaon-proton collisions are simulated in terms of multiple partonic soft and hard in- 
teractions, where each such partonic interaction leads to two QCD color strings, which subsequently 
fragment into hadrons. Diffraction dissociation is implemented as two-channel model of excited states 
for projectile and target particles, similar to the Good- Walker model flVl] of diffraction dissociation. 
The minijet cross section is calculated within the QCD-improved parton model using the parton den- 
sity parametrization of Gliick, Reya and Vogt II172II . Saturation is accounted for by introducing an 
energy-dependent transverse momentum cutoff for distinguishing soft and hard interactions II173I . 




Figure 10: Comparison of SIBYLL predictions with fixed target data from NA49 111 7411 at Eiah = 
158 GeV (left) and different data sets from colhder experiments (right) ^ = 54 - ISOOGeV 111701 . 

The scattering of hadrons off nuclei is modeled using the Glauber approximation. Realistic shell 
model distributions for the nucleon densities in the different nuclei are implemented 111751 . The semi- 
superposition model is used for the description of the interaction of nuclei as beam particles with 
air Iil76.l . The MC results are compared to fixed-target and collider data in Fig.fTOl 

Pomeron calculus: cross sections, diffraction and MC simulation 

Sergey Ostapchenko (NTNU, Trondheim) 

Monte Carlo (MC) generators of hadronic interactions are standard tools for data analysis in high 
energy collider and cosmic ray (CR) fields. A general framework for developing such generators is 
provided by the Reggeon Field Theory 111771 . The key technique is the one proposed by Abramovsky, 
Gribov and Kancheli (AGK) 1 178 1, which allows one to relate partial cross sections for various config- 
urations of hadronic final states to certain unitarity cuts of elastic scattering diagrams. The procedure 
generally requires the knowledge of the amplitude for an "elementary" elastic rescattering process 
("Pomeron exchange") which corresponds to an underlying microscopic parton cascade, of the ver- 
tices for Pomeron-hadron coupling, and of the vertices for Pomeron-Pomeron (PP) interactions, if 
non-linear interaction effects are to be taken into account. In principle, all the three ingredients may 
be specified using different approaches, ranging from purely phenomenological parametrizations to 
BFKL-based treatment. In practice, one usually considers "soft" and "semihard" contributions to the 
elastic rescattering process, depending on whether it is dominated by a purely nonperturbative soft 
parton cascade or parton evolution extended to moderately large virtualities While soft parton 
evolution is described by phenomenological parametrizations (e.g., as soft Pomeron emission), its 
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extension to higher \q^\ is typically treated within the DGLAP formalism. 

It is worth stressing that the starting point for a development of a MC generator is the derivation 
of the complete set of partial cross sections for various hadronic final states. The latter are defined 
by the number of elementary particle production processes ("cut Pomerons") and by a particular 
arrangement of those cut Pomerons in rapidity-impact parameter space. Thus, calculations of partial 
cross sections involve full resummation of virtual (elastic) rescatterings which are described as uncut 
Pomeron exchanges. Such an analysis becomes especially simple when using eikonal vertices for 
Pomeron-hadron coupling and neglecting Pomeron-Pomeron interactions, which leads to the standard 
eikonal description of the elastic scattering amplitude and to simple expressions for partial cross 
sections. In particular, one arrives to the Poisson distribution for a number of elementary production 
processes for a given impact parameter. 

However, the necessity to describe non-linear corrections to the interaction dynamics, related 
to parton shadowing and saturation, forces one to take Pomeron-Pomeron interactions into account, 
which significantly complicates the formalism. Indeed, with the energy increasing, so-called en- 
hanced {PP interaction) graphs of more and more complicated topologies start to contribute signif- 
icantly to the scattering amplitude and to partial cross sections for particular hadronic final states. 
Thus, dealing with enhanced diagrams, all order resummation of the corresponding contributions is 
a must, both for elastic scattering diagrams and for the cut diagrams representing particular inelas- 
tic processes. Secondly, it is far non-trivial to split the complete set of cut enhanced diagrams into 
separate classes characterized by positively-defined contributions which could be interpreted proba- 
bilistically and employed in a MC simulation procedure. While the first problem has been addressed 
in II 1791418^1 1, the MC implementation of the approach has been discussed in 11182] . Let us briefly list 
the main results of the analysis of Refs. II179I4I821 . 

• Non-linear contributions to the interaction dynamics aie not small corrections, but dominate 
the high-energy behavior of hadronic cross sections and particle production ||1791I180|[T82I . 

• Unlike inclusive particle (jet) cross sections, partial cross sections of hadronic final states which 
involve high transverse momentum jets can not be expressed via universal parton distribution 
functions (PDFs) measured in DIS experiments; rather they depend on "reaction-dependent 
PDFs" which involve parton rescattering on both the parent (e.g., projectile) hadron and the 
partner (here, target) hadron [180|. In other words, coUinear QCD factorization is inapplicable 
to exclusive hadronic final states. 

• An analysis of partial contributions of different classes of enhanced diagrams has shown that 
neither a resummation of contributions of "fan"-like graphs or of a more general class of "net"- 
like enhanced graphs nor of the ones of "Pomeron loop" diagrams alone is sufficient for a 
correct description of hadronic cross sections in the high-energy limit III8II . Instead, both 
classes of diagrams have to be taken into account. 

• Calculations of diffractive cross sections require a proper resummation of absorptive correc- 
tions to the contribution of diffractively cut sub-graphs, which seriously reduce the probability 
for a rapidity gap survival, in addition to the usual (eikonal) rapidity gap suppression (RGS) 
factor II181II . In particular, restricting oneself with the simplest triple-Pomeron contribution to 
single high mass diffraction cross section, one arrives to a contradiction with s-channel unitar- 
ity, even if the eikonal RGS factor is included. 

• The complete set of unitarity cuts of elastic scattering diagrams can be re-partitioned into 
a number of positively-defined contributions which define partial cross sections for certain 



32 



"macro-configurations" of tiie interaction (i^-cliannel unitarity) II182I . For eacli of tliose macro- 
configurations, the pattern of secondary particle production can be reconstructed in an iterative 
fashion using a MC procedure. 

One has to mention, however, that the approach of II 179141821 has a serious drawback of neglecting 
energy-momentum correlations between multiple scattering processes at the amplitude level |183l. 
Additionally, the discussed treatment uses phenomenological parametrization for Pomeron-Pomeron 
interaction vertices and neglects hard (high \q^\) PP coupling. Hence, the scheme is unable to de- 
scribe the dynamical evolution of the saturation scale in hadron-hadron scattering. 

Understanding the "ridge" in proton-proton scattering at 7 TeV 

Klaus Weme¥^ ( Subatech. Nantes) 

The CMS collaboration published recently results ||20| on two-particle correlations in Ar] and 
A(|), in PP scattering at 7 TeV. Most remarkable is the discovery of a ridge-like structure around 
Ar] = 0, extended over many units in Ar], referred to as "the ridge", in high multiplicity pp events. A 
similar structure has been observed in heavy ion collisions at RHIC, and there is little doubt that the 
phenomenon is related to the hydrodynamical evolution of matter. This "fluid dynamical behavior" is 
actually considered to be the major discovery at RHIC. 




Figure 11: Two-particle correlation function R versus Ar] and A(|) for high multiplicity events in pp 
collisions at 7 TeV, as obtained from a hydrodynamical evolution based on flux tube initial conditions. 
We consider particles with pj between 1 and 3 GeV/c. 

So does PP scattering provide as well a liquid, just ten times smaller than a heavy ion collision? It 
seems so ! We showed recently HI 841 that if we take exactly the same hydrodynamic approach which 
has been so successful for heavy ion collisions at RHIC fT85l, and apply it to pp scattering, we obtain 
already very encouraging results compared to pp data at 0.9 TeV. In this paper, we apply this fluid 
approach, always the same procedure, to understand the 7 TeV results. In Fig. [TT] we show that our 
hydrodynamic picture indeed leads to a near-side ridge, around A(|) = 0, extended over many units in 
Ar]. For the pure basic string model, without hydro evolution, one finds no ridge! This shows that the 
hydrodynamical evolution "makes" the effect. 

^^In collaboration with lu.Karpenko and T.Pierog. 
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It is easy to understand the origin of the ridge, in a hydrodynamical approach based on flux tube 
initial conditions, see [1861. Imagine many (say 20) flux tubes of small transverse size (radius ^ 0.2 
fm), but very long (many units of space-time rapidity r] , ). For a given event, their transverse positions 
are randomly distributed within the overlap area of the two protons. Even for zero impact parameter 
(which dominated for high multiplicity events), this randomness produces azimuthal asymmetries. 
The energy density obtained from the overlapping flux tubes shows an elliptical shape. And since the 
flux tubes are long, and only the transverse positions are random, we observe the same asymmetry at 
different longitudinal positions. So we observe a translational invariant azimuthal asymmetry! 

If one takes this asymmetric but translational invariant energy density as initial condition for 
a hydrodynamical evolution, the translational invariance is conserved, and in particular translated 
into other quantities, like the flow. At a later time, at different space-time rapidities, the flow is 
more developed along the direction perpendicular to the principal axis of the initial energy density 
ellipse. This is a typical fluid dynamical phenomenon, referred to as elliptic flow. Important for this 
discussion: the asymmetry of the flow is again translational invariant, the same for different values 
of r\s. Finally, particles are produced from the flowing liquid, with a preference in the direction of 
large flow. This preferred direction is therefore the same at different values of r],. And since r\s and 
pseudorapidity r\ are highly correlated, one observes a Ar\A(^ correlation, around A(|) = 0, extended 
over many units in Arj: a particle emitted a some pseudorapidity r\ has a large chance to see a second 
particle at any pseudorapidity to be emitted in the same azimuthal direction. 

FLUKA Monte Carlo 

Maria Vittoria Garz.ell^^ (INFN Milano & Univ. Granada) 

The FLUKA Monte Carlo code 111871 . also interfaced with the dpmjet code II188II for the treatment 
of nucleus-nucleus interactions, is being used in cosmic-ray physics. More detailed information on 
the physics models for h-h, h-A, A— A interactions relevant for cosmic-ray Physics adopted in FLUKA 
and DPMJET is available in the literature |189]. Experimental observables involving muons, like 
cosmic /^^/^^ charge ratios and the ^ decoherence function of bundles detected underground, can 
be used to test the hadronic interaction models. 

Cosmic liT charge ratios have been measured both by detectors at accelerator sites (L3-I-C, 
CMS, preliminary data from ALICE) and by passive underground experiments (Utah, MINOS, OPERA). 
The results of FLUKA simulations on ^^/^^ charge ratio are shown together with available experi- 
mental data 111 901 in Fig. [12] (left). They turn out to be completely compatible with the CMS and 
L3-I-C data, whereas they slightly underestimate the MINOS data (however not completely confirmed 
by the OPERA ones, especially at the highest energies). A possible reason of this discrepancy is the 
fact that the /K^ charge ratio or the ratio between /K^ and 7i+/7l^ charge ratios in FLUKA can 
be underestimated. The differences in the steepness of xpeynman distributions of and K is also 
fundamental. To test these hypotheses new data on K and 71 production at high-energy accelerators 
are urgently needed. 

The decoherence function, i.e. the distribution of the average distance between ^ pairs in a bundle 
detected underground, is another observable that allows one to test the transverse structure predicted 
by the hadronic interaction models used to interpret the results of experiments on cosmic-rays. Such 
a measurement is highly sensitive to the pT distributions of % and K produced in the first stages of the 
shower formation process. We have verified that the decoherence function is robust against a change 
of the details of the cosmic -ray primary spectrum. FLUKA -i- DPMJET reproduces the experimental 

^^G. Battistoni, M.V. Gai'zelli, A. Margiotta, S. Muraro, M. Sioli for the FLUKA Collaboration 
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data from the 



MACRO experiment 111911 . as shown in Fig. [T2l (right). 




e, cos 0" (GeV) Separation (m) 

Figure 12: Left: fj'^ /fJ^ charge ratio as a function of the surface energy E^cosQ: FLUKA simulation 
(pink bands representing an average charge ratio over all E^cosQ values in the corresponding abscissa 
interval) vs. experimental data (Utah, LVD, L3+C, MINOS, CMS and OPERA). Right: decoherence 
distribution of /j pairs in a bundle (histogram: FLUKA+DPMJET; symbols: MACRO data) 

Comparison of model predictions to LHC data 

Tanguy Pierog (KIT, Karlsruhe) 

In April 2010, the ALICE HI 921 collaboration published for the first time the pseudorapidity 
distribution of charged particles at a center-of-mass energy of 7 TeV [26|. The analysis is based on 
a trigger called lnel>0 which consist of having at least one particle with |ri| < 1. Unlike the Non 
Single Diffractive trigger, this does not include any model-based correction and allows then for an 
easy comparison with any hadronic interaction models. In Fig. [13] we compare the ALICE data at 
900 GeV, 2.36 and 7 TeV with the models commonly used in air shower simulations, namely QGSJET 
01 HH and n-03 HH, SIBYLL 2.1 UTOJ and epos 1.99 [195J, using the same trigger. 
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7 - QGSJETII 7000 GeV 
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Figure 13: Pseudorapidity distribution of charged particles measured by ALICE m pp {Inel > 0) 
colhsions compared to the predictions of QGSJET 01 and 11-03, SIBYLL 2.1, and EPOS 1.99 models. 

Even if none of these MCs can reproduce perfectly all the data, it is interesting to notice that 
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the experimental results are well bracketed by the models used for cosmic ray analysis. The two 
extremes are EPOS 1.99 which underestimate the growth of the particle density and QGSJET 11-03 
which overestimates it. As a consequence we can say that the predictions of the Monte-Carlo models 
based on Regge field theory used for extended air shower analysis are compatible with the recent 
data from LHC, and hence do not show any tendency of dramatic change in minimum bias hadronic 
physics. Detailed results and dependence on trigger conditions not shown here indicate that soft 
physics is not yet completely understood. 
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